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ABSTRACT 

Aims. We study possible signs of asymmetry in the bright optical transient in NGC 300, to obtain independent information about the 
explosion mechanism, the progenitor star and its circumstellar environment. 

Methods. Using VLT-FORS1, we obtained low-resolution optical linear spectropolarimetry of NGC 300 OT2008-1 on two epochs, 
48 and 55 days after the discovery, covering the spectral range 3600-9330A. 

Results. The data exhibit a continuum polarization at a very significant level. At least two separate components are identified. The 
first is characterized by both a strong wavelength dependency and a constant position angle (68.6±0.3 degrees), which is parallel to 
the local spiral arm of the host galaxy. The latter is aligned along a completely different position angle (151.3±0.4). While the former 
is identified as arising in the interstellar dust associated with NGC 300, the latter is most likely caused by continuum polarization by 
dust scattering in the circumstellar environment. No line depolarization is detected in correspondence to the most intense emission 
lines, disfavoring electron scattering as the source of intrinsic polarization. This implies that there is a very small deviation from 
symmetry in the continuum-forming region. Given the observed level of intrinsic polarization, the transient must be surrounded by a 
significant amount of dust (>4xl0~ 5 M Q ), asymmetrically distributed within a few thousand AU. This probably implies that one or 
more asymmetric outflow episodes took place during the past history of the progenitor. 
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1. Introduction 

Many surveys have focused their interest on poorly known ob- 
jects lying in the gap in the absolute magnitude distribution that 
separates luminous novae (Mr > -10) and faint supernovae 
(M R < -14) (e.g., Rau et al. 120081 ). Although this gap had been 
understood to be populated mostly by eruptions of luminous 
blue variables (LBVs, see Maund et al. [2006; Van Dyk et al. 
2006), it is now evident that new types of explosions may be 
responsible for the transient events occurring in this luminos- 
ity range, from exotic (LBV-like) outbursts of Wolf-Rayet stars 
(Pastorello et al. l2007at to other long duration transients, such as 
M85 OT2006-1 (Kulk arni et a l.[200ll Rau et al. l2008l Ofek et al. 
2008 ; Pastorello et al. 2007b), whose nature remains debated. In 
this context, a significant impulse to the study of objects in this 
gap was brought about by the discoveries of two nearby objects, 
namely SN 2008S in NGC 6946 and the 2008 luminous optical 
transient in NGC 300 (hereafter NGC 300 OT2008-1). 

NGC 300 OT2008-1 was discovered on May 14, 2008 
(Monard 2008). A spectrum obtained the day after, revealing 
Ha, H/3, Ca II near-IR triplet, forbidden [Ca II] doublet in emis- 
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* Based on observations made with ESO Telescopes at Paranal 
Observatory under Program ID 281.D-5016. 



sion, and Ca II H&K absorptions (Bond, Walter & Velasquez 
120081) , shared striking similarities with that of SN 2008S (Prieto 
2008 ). The photometric and spectroscopic properties of the tran- 
sient in NGC 300 led to the conclusion that this was not a classi- 
cal nova, an LBV, or a supernova (Bond et al. 2009). According 
to Bond et al. (2008 , 2009 ), the spectrum was reminiscent of that 
of V838 Mon during the early phases of its outburst (Wisniewski 
et al. 120031 ). Broad-band optical and near-IR photometric, plus 
low and intermediate resolution spectroscopic follow-up obser- 
vations were presented by Bond et al. (2009 ), while Berger et al. 
(2009 ) reported on UV, radio, and X-ray imaging and high reso- 
lution spectroscopy. Remarkably, no source was detected at the 
position of NGC 300 OT2008-1 in the radio and X-ray observa- 
tions, which indicates that the transient was at least one order of 
magnitude (in the radio domain) and two orders of magnitudes 
(in the X-rays region) fainter than any SN discovered to date. In 
addition, high-resolution spectroscopy suggests the presence of 
a complex circumstellar environment resulting either from pre- 
vious asymmetric ejections from the progenitor or from a com- 
panion wind (Bond et al. l2009t Berger et al. l2009l . 

An inspection of pre-explosion, archival HST images re- 
vealed that no star was detectable at the outburst position plac- 
ing stringent upper limits on the progenitor's luminosity, which 
was interpreted as an additional argument in favor of a non- 
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supernova origin (Berger & Soderberg 12008b . However, a pro- 
genitor candidate was detected in pre-outburst Spitzer images 
(Prieto 2008 ; Berger et al. 2009 ), and it was identified as a mas- 
sive star (15-20 M ), with a spectral energy distribution very 
similar to that exhibited by the progenitor of SN 2008S (Prieto 
et al. 120081 Botticella et al. 2009). By means of an independent 
method based on the study of the surrounding stellar population, 
a similar range for the progenitor mass (12-17 M ) was found by 
Gogarten et al. (120091) . while Bond et al. ([2009 ) gave a slightly 
lower mass range (10-15 M ). An analysis of the Spitzer/IRS 
mid-IR spectrum of NGC 300 OT2 008-1 o btained three months 
after the outburst led Prieto et al. (120091 ) to conclude that, al- 
though the presence of a massive star cannot be definitely ex- 
cluded, a lower mass (6-10 M , carbon-rich AGB/super-AGB, 
or a post-AGB) precursor would provide a better match to their 
mid-IR observations. 

Albeit with some scatter in the estimate of the progenitor's 
mass, most studies favor an exotic eruption of a moderate-to- 
massive star for NGC 30 OT2008-1 (and SN 2008S). However, 
Botticella et al. (12009) emphasized the overall similarity be- 
tween the light curves of these transients and those of type II-L 
SNe. The main argument was the late-time flattening, which is 
consistent with the slope expected from the radioactive decay of 
56 Co into 56 Fe. If this mechanism plays a role in powering the 
light curve of these transients, then a SN explosion is the most 
likely engine (electron-capture SNe, see Botticella et al. 2009; 
Pumo et al. 2009, submitted). 

Soon after the discovery, we started an optical/near-IR 
follow-up using a number of ground-based facilities. The re- 
sults will be presented and discussed in a forthcoming paper 
(Pastorello et al. 2009, in preparation). In this paper we focus 
on the spectropolarimetry of NGC 300 OT2008-1 that we ob- 
tained in July 2008 with VLT-FORS 1 on two epochs separated 
by a week, with the aim of detecting possible signs of asymme- 
tries in the ejecta and/or in the circumstellar environment that 
could provide additional constraints on the progenitor's nature. 

The paper is organized as follows. In Sect. |2]we discuss the 
observations and the data reduction techniques. In Sect. [3] we 
present the flux spectra, while Sect. |4] deals with the spectropo- 
larimetric data sets. In Sect. [5] we discuss the results, and we 
summarize our conclusions in Sect. [6] Appendix |A| contains an 
analysis of the effects of second order contamination on spec- 
tropolarimetry and its application to the data presented in this 
work, while Appendix [B] describes the effect of multiple weak 
polarizers on the resulting polarization signal. 



2. Observations and data reduction 

We observed the OT in NGC 300 on 2 di fferent epochs, 48 
and 55 days after its discovery (Monard 12008b . using the 
FOcal Reducer/low-dispersion Spectrograph (hereafter FORS1), 
mounted at the Cassegrain focus of the ESO-Kueyen 8.2m tele- 
scope (Appenzeller et al. 11998b . In this multi-mode instrument, 
equipped with a mosaic of two 2048x4096 pixel (px) E2V 
CCDs, polarimetry is performed by introducing into the optical 
path a Wollaston prism (19" throw) and a super-achromatic half- 
wave plate (HWP). To reduce some known instrumental prob- 
lems (see Patat & Romaniello 2006a), we always used 4 half- 
wave plate angles (#,-=0, 22.5, 45 and 67.5 degrees). All spec- 
tra were obtained with the low-resolution G300V grism coupled 
with a 1 .0 arcsec slit, giving a spectral range 3200-9330 A, a dis- 
persion of ~3.2 A px -1 , and a resolution of 1 1.6 A (FWHM) at 
5800 A. To cover the blue part of the optical spectrum (<4350A), 
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Fig. 1. Flux spectra of NGC 300 OT2008-1 on 01-07-2008 (up- 
per) and 08-07-2008 (lower). For presentation, the spectrum of 
the second epoch was shifted by -0.4. The main spectral features 
are identified. 

where significant Ca II H&K absorption is clearly evident, we 
have not used an order sorting filter. Since the object is rather 
red (B - V ~1.1, see Sect. [3}, the second order contamination, 
which for this grism starts at about 6600 A (O'Brien 12008b . is 
expected to be small and its effects on the spectropolarimetry of 
NGC 300 OT2008-1 negligible. This is discussed more quanti- 
tatively in Appendix [A] 

Data were bias and flat-field corrected, and wavelength cali- 
brated using standard tasks within IRAfQ The RMS error in the 
wavelength calibration is about 0.7A. The ordinary (upper) and 
extraordinary (lower) beams were processed separately. Stokes 
parameters, linear polarization degree, and position angle were 
computed by means of specific routines written by us. Finally, 
polarization bias correction and error estimates were performed 
following the prescriptions described by Patat & Romaniello 
(2006a), while the HWP zeropoint angle chromatism was cor- 
rected using tabulated data (O'Brien 2008). To increase the 
signal-to-noise ratio, the final Stokes parameters were binned in 
~25.8A wide bins (8 pixels). This turns into an RMS error in the 
polarization of 0.1% at 6000A, where the continuum signal-to- 
noise ratio per resolution element reaches its maximum (720 and 
530 on the first and the second epoch, respectively). The spec- 
tra are severely affected by fringing above ~8000A, so that the 
observed noise at the red edge of the wavelength range is sig- 
nificantly greater than the statistical formal errors. The instru- 
ment stability (including evolution of the fringing pattern) was 
checked by comparing £q, e\j pairs (Maund 12008b at the two 
epochs. Their average differences 

angleAe) = (eg - eu) were found to be 0.07+ 0.02% and 
0.13±0.02% on the first and the second epoch, respectively. 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, under contract with the National Science 
Foundation. 
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Fig. 2. Observed spectropolarimetry of NGC 300 OT2008-1 at the two epochs on the Stokes Q-U plane. For presentation, the data 
were binned to 16 pixels (57A). The dotted circles are placed at constant polarization levels (0.5, 1.0, 1.5, and 2.0%). The dashed 
line is a weighted least squares fit to the data in the displayed wavelength range. For comparison, the colored line sequence of dots 
in the upper left quadrant traces a Serkowski law (Serkowski, Matheson & Ford 1975) for /i max =4800A, P max =lA%, and /T=0.81 
with an arbitrary position angle 0=60 degrees. The circled dot indicates the position of (Q\, U\) corresponding to the minimum 
polarization solution (see Sect. |4. 1]>, while the vector indicates the allowed region for (gi, U\). 



Table 1. Log of the VLT-FORS1 observations of 
NGC 300OT2008-1. 



Date 


MJD 


Exp. Time 


Airmass 


(UT) 




(s) 




01-07-2008 


54648.4 


2x500+2x400 


1.1 


08-07-2008 


54655.3 


660+3x500 


1.2 



These values are within the measurement errors per resolution 
element (typically <0.1%), thus confirming the stability of the 
instrument in the relevant time interval. However, they also im- 
ply possible systematic errors of up to about 0. 1 %, which can be 
considered as the maximum accuracy one can reach with FORS 1 
using 4 HWP positions (see also Patat & Romaniello 2006a). 

Flux calibration was achieved by observing a spectrophoto- 
metric standard star and inserting full polarimetric optics (HWP 
angle set to 0). A log of the observations is given in Table[T] 

3. Flux spectra 

N GC 300 OT2008-1 was fo llowed up in great detail by Bond et 
al. (J2009 ) and Berger et al. ((2009 ), to which we refer the reader 
for a comprehensive account of the spectroscopic properties of 
this object and its evolution. In the following, we provide a brief 
description of the flux spectra, while they will be discussed in 
more detail in the analysis of the whole spectroscopic data set of 
NGC 300 OT2008-1 (Pastorello et al., in preparation). 

The spectra, obtained by coadding for each epoch the ordi- 
nary and extraordinary beams for all the four HWP angles, are 



presented in Fig. [T] Between the two epochs, which are sepa- 
rated by seven days, no evolution is evident and the spectrum 
consists of a red continuum with a number of superimposed 
narrow emission lines. As Bond et al. (2008) pointed out, the 
most prominent spectral lines can be identified with the Balmer 
series of Hydrogen, Ca II 8498A, 8542A, 8662A, and [Ca II] 
7291, 7324A. The H lines can be identified to wavelength as 
sort as H6 and the Balmer jump at ~3650A is clearly detected 
(Fig. [I}. Additionally, our spectra appear to contain He I 3888 A, 
4571A, 5015A, 5876A, 6678A, 7065A absorption lines, which 
are clearly seen in the higher resolution spectra presented and 
discussed by Berger et al. ( 120091 . Finally, Ca II H&K and Na ID 
are visible in absorption and are contaminated by unresolved 
features arising in the ISM (see also Berger et al. 2009). 

While the absolute flux scales should be compared to broad- 
band photometry when this becomes available, the overall shape 
of the spectra is expected to be accurate to within 10-15%. 
The colors deduced from synthetic photometry are B - V=+l.l, 
V - R-+0.6, and V - 7=+ 1.1 for both epochs, with RMS un- 
certainties of about 0.15 mag. These values are fully compatible 
(within the quoted errors) with the broadband photometry re- 
ported by Berger et al. d2U09l for June 19 (B - V=1.01 + 0.03; 
see their Table 2). Because of the relatively high value of B - V, 
the flux contamination by the second order above 6500A is ex- 
pected to be less than a few percent (see Appendix [A|. 

The measured FWHMs of Ha, UB, and Ca II 8662A are 
18.6+0.3A, 13.8+0.2A and 23.8+0.4A, respectively. Once cor- 
rected for the instrumental resolution (11. 6 A), these values cor- 
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Fig. 3. Differences in Q (upper panel) and U (lower panel) 
Stokes parameters between the two epochs (second minus first). 
Empty symbols indicate the data points rejected by the k-cr clip- 
ping. 

respond to expansion velocities of 660, 460, and 720 km s _1 . In 
general, the spectrum of NGC 300 OT2008-1 is very similar to 
those reported for SN 2008S by Smith et al. J2008b . Based on 
the very low luminosity and expansion velocities, these authors 
concluded that the object previously designated as SN 2008S 
is actually a SN impostor (Van Dyk et al. 2000), similar to the 
giant eruption of a luminous blue variable, and shares the spec- 
tral properties of the Galactic hypergiant IRC+ 10420 (Jones et 
al. 119931 1. Here we note that, in general, SN 2008S differs from 
a classical impostor, in that SN 2008S has a linear and slow 
evolving luminosity decline, a spectrum that matches those of 
LBV eruptions well, excluding its very strong [Ca II] doublet. 
Based on similar considerations, both Bond et al. d20091 > and 
Berger et al. ( 120091 concluded that NGC 300 OT2008-1 was 
the eruption of a dust-enshrouded, massive star. An alternative 
scenario was proposed by Prieto et al. ( 120081 1, Thompson et al. 
d2009l , and Botticella et al. ( 120091 , who argued that SN 2008S 
and NGC 300 OT2008-1 were generated by the explosion of an 
electron-capture SN. 



Fig. 4. A DSS image of NGC 300. The location of the optical 
transient is marked by a cross. The arrow indicates the polar- 
ization position angle of the dominant axis (the average of the 
values measured during the two epochs was used: 0=68.6 de- 
grees). 

nant axis. This does not pass through the origin of the Q - U 
plane, implying that there must be an additional component with 
a polarization wavelength dependence which differs from that 
along the dominant axis (see Appendix IB). We also note that 
the two data sets span similar ranges of values along the domi- 
nant direction, but that a rigid shift is present between the two 
epochs. This is more clearly illustrated in Fig. [3] where we 
plot the differences AQ and AU between the two epochs. The 
clipped average differences (computed in the wavelength range 
3900-9300A) are <A0=-O.l 1+0.02% (191 spectral bins) and 
(AU)=+0. 14+0.02% (193 spectral bins), and the RMS deviation 
from the average is 0.23 and 0.21% for the two Stokes parame- 
ters, respectively. Given that the observed grey variation is very 
small and comparable to the instrument stability between the two 
epochs (~0. 1%, see Sect. [2]), it is unclear whether this is real or 
caused by an incompletely removed instrumental effect. 

To characterize the dominant axis, we fitted the data with a 
linear relation 



4. Spectropolarimetry 

The spectropolarimetric data for NGC 300 OT2008-1 are pre- 
sented in the Q - U plane in Fig. [2] The object exhibits a sig- 
nificant polarization, which reaches a maximum level of about 
1.4%. During both epochs, the data cover a range of polariza- 
tion values which are definitely much higher than expected for 
a typical Serkowski law (Serkowski et al. 1975; see also next 
section) with a similar maximum polarization (see Fig.|2]and its 
caption). No obvious change in polarization is seen across the 
most prominent emission lines (Sect. 4.2 1. 



U = U {) + Q tan26» 



(1) 



4. 1 . Continuum polarization 

The most remarkable aspect of the Q - U plot is the align- 
ment of the data points along a straight line, which is indica- 
tive of a polarization component with a clearly defined domi- 



The results of the least squares fitting for the two epochs in the 
wavelength range 4000-9000A are reported in Table [2] (upper 
part), which includes the RMS deviation from the best-fit rela- 
tion (cr), the Pearsons linear correlation coefficient (r ry ) and x 1 
per degrees of freedom (DOF). To exclude the most deviant data 
points, we performed a preliminary fit to all data, applied a 1.7- 
sigma rejection and finally refit the remaining data. The large 
correlation factors confirm the presence of a dominant axis at a 
high significance level. This indicates that the main source of po- 
larization is interstellar extinction. Interstellar polarization (ISP), 
caused by the presence of aligned asymmetric dust grains along 
the line of sight, is characterized by a wavelength-independent 
position angle (Serkowski et al. 1 1975 1 Whittet et al. 1992). This 
conclusion is supported by the observed dominant axis being 
tangential to the spiral arm of NGC 300 close to the outburst 
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Fig. 5. Minimum polarization solution for 08-07-2008. The 
dashed line traces the observed dominant axis, while the point 
(Qi< U\) indicates the wavelength-independent component. The 
empty circles mark the rejected data points. 



Table 2. Dominant axis least squares fit parameters. 



Date 


Uo 


e 


(X 


r„ 


X 2 /DOF 




(%) 


(deg) 


(%) 






01-07-2008 


-0.42 (0.02) 


69.6 (0.4) 


0.19 


-0.94 


432/180 


08-07-2008 


-0.46 (0.02) 


67.6 (0.5) 


0.23 


-0.92 


333/179 


01-07-2008 


-0.43 (0.02) 


68.6 (0.4) 


0.20 


-0.94 


463/180 


08-07-2008 


-0.45 (0.02) 


68.6 (0.5) 


0.22 


-0.92 


311/179 



Note: Values in parenthesis are the statistical RMS uncertainties. 
The lower part of the table reports the same values after rotation 
to the average dominant axis (0=68.6 degrees). 



position, as shown in Fig. |4] Since dust grains tend to be aligned 
in the direction of the galactic magnetic field, which in turn fol- 
lows the spiral pattern (see Scarrot, Ward-Thompson & Warren- 
Smith 1987), this strongly indicates that not only the observed 
polarization is caused by directional extinction, but also that the 
polarization arises within interstellar material associated with a 
spiral arm of NGC 300. Similar alignments were detected using 
spectropolarimetry in two supernovae, namely SN 200 lei (Wang 
et al. EUOll ) and SN 2006X (Patat et al. 120091 . From now on we 
therefore assume that this component is constant in time. 

The estimated position angles at the two epochs differ by 
2.0+0.5 degrees. Since this variation is only marginally signifi- 
cant, we adopt the average of the two values (68.6±0.5 degrees) 
as the representative value. For consistency, we rotated the two 
data sets by -1.0 and +1.0 degrees, respectively. The results of 
the least squares fitting after this operation are shown in Table [2] 
(lower part). The recomputed average differences in Q and U are 
<A0=-0. 12+0.02% and <AJ/>=+0.1 1+0.02%. This means that 
the systematic, wavelength-independent shift observed between 
the two epochs is not related to this small rotation. 

As we said, the position of the points on the Q - U plane 
requires an additional polarization component. The simplest so- 
lution is given by a wavelength-independent polarization, repre- 
sented by one single point (Q\ , U\ ) on the Q - U plane. In these 
circumstances, the total, observed Stokes parameters Q and U 




4000 6000 8000 

Observed Wavelength (A) 

Fig. 6. Dominant (upper panel) and orthogonal (lower panel) 
polarization on 01-07-2008 after subtracting the wavelength- 
independent component Pi. For presentation, the data were 
binned to 16 pixel (57A). The solid curve in the upper panel 
is a Serkowski law with P mflA =4.8%, /l mfl . v =4250A, and K-3.0. 



are given by Q = gi + Qi{X) and U = U\ + U%(?i), where Q 2 
and Ui are the Stokes parameters that characterize the polariza- 
tion aligned along the dominant axi^] Of course, if the ratio of 
Q2 to U2 must remain constant (as required by the wavelength- 
independent position angle), Q\ and U\ must obey Eq. [T| i.e. 
(fii> U\) must be aligned along the observed dominant axis. One 
additional constraint related to the definitions of Q and U is that 
all the (Q2, U2) points have to be confined to one single quadrant 
of the Q—U plane (Wang et al. 120031 . which, in our case, is that 
with Q < and U > 0. This imposes that Q\ is higher than a 
certain value Q„,i„, but does not prevent Q\ being any higher ar- 
bitrary value (see the arrow in Fig. [2]). For this reason, although 
constrained, the solution is not unique. As a consequence, while 
the position angle of the dominant axis (#2) is fixed, the corre- 
sponding polarization P2(/l) - ^Q\(X) + i/|(/l) is undefined up 
to an additive constant, which depends on the particular choice 
of the wavelength-independent polarization Pi = + U\- 

Since (Q\ , U\) must lie along the dominant axis, the correspond- 
ing value of 9\ = arctan U\IQ\ is allowed to vary only within a 
rather limited range (smaller than 10 degrees), regardless how 
large the value of Q\ is. 

4.1.1. Minimum polarization solution 

In the following, we consider what we will indicate as minimum 
polarization solution (MPS), i.e., the solution one obtains for 
Q\ = Qmim an d that leads to the minimum possible values of the 
polarization degree for both components. Given the wavelength 
dependency shown by the data (Fig. [2}, this means that the MPS 
corresponds to the case where the polarization is close to zero at 
the red edge of the observed wavelength range (~9300A). The 



2 In the weak polarization approximation the components are vecto- 
rially additive. See AppendixlB] 
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Table 3. Minimum polarization solutions for the two epochs. 



Date 


Gi 




Pi 


0i 




(%) 


(%) 


(%) 


(deg) 


01-07-2008 


+0.75 


-1.12(0.03) 


1.35(0.03) 


151.9(0.5) 


08-07-2008 


+0.63 


-1.03(0.03) 


1.21(0.03) 


150.7(0.5) 



Values in parenthesis are the statistical RMS uncertainties. 



solution was found as the value of Q\ that produces the min- 
imum constant value of P2W at any given wavelength, while 
fulfilling the constraints described above. Since we had shown 
that the two data-sets differ by a wavelength-independent quan- 
tity (cf. Fig. [3J, we computed gi f° r the second epoch and 
set the value for the first epoch such that A(2i=-0.12% (see 
above). Within the errors, this ensures that P2 W does not change 
between the two dates. The results of this procedure are pre- 
sented in Tableland illustrated in Fig.[5]for the second epoch. 
What emerges from this analysis is that the additive, wavelength- 
independent component is polarized at a highly significant level, 
which reaches at least 1.2%. Its position angle (151.3+0.4 de- 
grees, averaged over the two epochs) is very different from that 
of the dominant component, indicating that they must arise from 
dissimilar sources (see the discussion in Sect. [5). 

When, as in this case, there is a marked dominant direction, 
it is useful to compute the components parallel (Pj) and orthog- 
onal (P ) to this axis, as suggested by Wang et al. (2003 ). This 
is achieved by the following roto-translation in the Q— U plane: 

Pd = (Q- Qi)cos26 z + (U -Ui)sin20 2 , 
P„ = (C/-C/i)cos26» 2 -(e-ei)sin20 2 . 

While Pd is undefined up to an additive constant (dependent on 
the particular choice of Q\), P is not. Indeed, P a is purely re- 
lated to the deviations from the dominant axis. The decomposi- 
tion is shown in Figs. [6] and [7] for the two epochs. Two important 
results emerge from these plots: i) Pj has a smooth wavelength 
dependence, reaching a maximum (2.6±0.1%) at about 4900A 
and decreasing on either sides of the peak much faster than a 
Galactic Serkowski law (Serkowski et al. 11975b : ii) P is con- 
sistent with a null value. The RMS deviations in the wavelength 
range 4000-9000A are 0.13% and 0.16% during the two epochs, 
respectively. These values are fully compatible with the esti- 
mated measurement errors, and thus the two components Pi and 
P2 are sufficient to explain all the observed variance in the resid- 
uals. During the second epoch, there might be an indication of a 
mild, smooth wavelength dependency in P a (see Fig. [7] lower 
panel), but the statistical significance appears to be marginal. 
The same is true for a small bump seen at about 7800A during 
the first epoch (Fig. [6] lower panel). 

For highly reddened Galactic stars, Serkowski et al. (11975b 
showed that the linear polarization degree changes with wave- 
length according to the empirical law 

P(X) = P max exp(-tf In 2 ^p), 

where A max is the wavelength that corresponds to the maximum 
polarization P max - Witthet et al. ( 1 992 ) also demonstrated that 
K = 0.01 + 1.66 A max , where A max is expressed in microns. 
Typical values of K range from 0.7 to 1.4. A comparison be- 
tween the NGC 300 OT2008-1 data and the behavior expected 
on the basis of the Serkowski law immediately shows that the ob- 
served wavelength dependence is much steeper (Fig.|2|, indica- 
tive of values of K higher than 5 (see Figs.|6]|7]i. For comparison, 
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Fig. 7. Same as Fig.^for 08-07-2008. 

Winsniewski et al. (2003) report K ^1 for the Galactic, pecu- 
liar variable V838 Monocerotis, for which they could reproduce 
the observational data with a Serkowski law (see their Fig. 5). 
Interestingly, Pd is reasonably well fitted by a Serkowski law 
with P max =2.6+Q.l%, ,W=4880+40A, and #=5.6+0.1 (the 
Witthet et al. relation gives #=0.82 for the observed A max ). 

This marked deviation indicates that the dust mixture re- 
sponsible for the polarization most likely differs from a typi- 
cal Galactic mixture. Therefore, Serkowski's Galactic relation 
between color excess and maximum polarization (P max (%) > 
9 E B -v) probably does not hold for NGC 300 OT2008-1. In this 
respect, we note that the ISP wavelength dependency observed 
in several supernovae also clearly deviates from the Galactic 
Serkowski law (Leonard & Filippenko 120011 Leonard et al. 
120021 Maund et al. 120071 Patat et al . 120091 . For the hypothesis 
that the dominant polarization detected in NGC 300 OT2008- 
1 is produced by dust extinction, as all elements seem to indi- 
cate, one can attempt to estimate reddening from the observed 
level of polarization. Blindly using P max -2.6% in the quoted re- 
lation, one obtains an upper limit to the color excess given by 
Eb-v <0.3. The observed deviation from a Galactic Serkowski 
law casts serious doubts on the reliability of this estimate, and 
what we quote here is just for the sake of completeness. To our 
knowledge no polarimetric study on NGC 300 exists in the lit- 
erature. Therefore, it is not possible to tell whether the observed 
wavelength dependency is typical of the host galaxy or is a pe- 
culiarity of the line of sight to NGC 300 OT2008-1. However, 
we note that the extinction law along the line of sight to young 
clusters in NGC 300 was reported to be highly variable (Roussel 
etal.[2005]>. 

Based on photometric considerations, Bond et al. 
(120091 report £' B _y=0.4, while Berger et al. ( [2009] ) derive 
£^=0.05+0.05 from the equivalent width of the Na I D 
lines. They conclude that this low value for the reddening is 
consistent with the lack of obvious evidence of extinction in 
their optical and UV spectra, and point out that this is in line 
with evidence of dust destruction in the immediate surroundings 
of the outbursting star. The extinction arising within the Milky 
Way is very small (£Vy=0.013, Schlegel, Finkbeiner & Davis 
119981 ) and the bulk of reddening has to occur within the host 
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Table 4. Synthetic broadband polarimetry of NGC 300 OT2008- 
1. 



Epoch 


B 




V 




R 






I 




P(%) 


e 


P(%) 


e 


P(%) 


e 


P(%) 


9 


#1 


1.21 


79.0 


1.04 


77.5 


0.55 


84.6 


0.67 


141.7 


#2 


1.36 


76.4 


1.18 


74.1 


0.71 


78.2 


0.53 


133.8 


#2(*) 


1.21 


78.4 


1.02 


76.1 


0.55 


82.3 


0.67 


139.4 



(*) Recomputed after applying < AQ >=-0.11,< AU >=+0.14. 
Note: The RMS errors on polarization levels and position angles are 
0.02% and 0.3 degrees, respectively. 



galaxy, and so is most likely the case for interstellar polarization. 
The compilation of Heiles (1999 ) contains a few stars close 
in galactic coor dinates to NGC 3OC0 All of them have linear 
polarizations <0.1%, in close agreement with the low Galactic 
extinction along this line of sight. 

The best-fit global solution is compared to the observed 
Stokes parameters and position angle in Fig. [8 for the data ac- 
quired on 2008 July 8. The fit is fairly good for this epoch 
(y 2 =l.l, RMS=0.20%), while for the first epoch the observed 
data points show higher scatter redwards of ~7000A (^ 2 =2.8, 
RMS=0.28%). 

The synthetic B, V, R, and / broad-band polarizations are re- 
ported in Table [4] As we have shown for the spectra, the ob- 
served variations (all highly significant) can be explained in 
terms of evolution in the wavelength-independent component. 
After applying the systematic differences (AQ) and (AC/) to the 
second epoch and re-computing the synthetic broad-band po- 
larimetry, both the polarization angle and degree are consistent 
(within the errors) with those derived for the first epoch (Table|4] 
last row). Interestingly, the variation in Q\ and U\ takes place 
along a direction (PA=154. 1+3.2) consistent with the average 
position angle of P\ (151.3+0.4). If this variation were real, 
this would imply that the wavelength-independent component 
evolves just within its polarization level, while the position angle 
remains practically unchanged. This change in the polarization 
level is indicative of an intrinsic origin. However, since the ob- 
served variation is within the typical instrumental stability limits 
(~0.1%, see Sect.|2]i, no solid conclusion about its nature can be 
reached. 



4.2. Line polarization 

Although the resolution of our spectra is rather low (~12A 
FWHM), we inspected the data searching for possible line ef- 
fects corresponding to the most prominent emission lines. In do- 
ing this, we first vectorially subtracted the contribution of the 
component Pi (which we assumed to be of interstellar origin 
and to be constant across the line profile) directly estimated fit- 
ting the adjacent Qz(A) an d Uz(A) data points and interpolat- 
ing to the line wavelength^] The case of Ha, the most intense 
feature (see Fig. [TJ, is presented in Fig. [9] for the first epoch. 
Across the line profile there is no evidence of the polarization 
changes that were previously detected, for instance, in Nova 
Cygni 1992 (Bjorkman et al. [T994l) . Nova Scuti 2 003 ( Kawabata 
et al. l2006l l. and V838 Mon (Wisniewski et al. 120031 . A sim- 
ilar behavior is observed during the second epoch and for all 
other prominent emission lines (HB and Ca II NIR triplet), in the 



sense that the observed polarization remains at the wavelength- 
independent level of Pi. If the emission line is essentially unpo- 
larized (Harrington & Collins 119681 ). the ratio of line to contin- 
uum polarization is given by the following semi-empirical rela- 
tion (McLean & Clarke H979]): 



P(A) 



1 



1 + F(A)/F C 



3 http://vizier.cfa.harvard.edu/viz-bin/VizieR?-source=II/226 

4 Note that this is a pure shift of Q and U that does not modify their 



shape. 



where F(A) is the line intensity profile and F c is the adjacent 
continuum level. Based on this relation and given that, at the 
peak of Ha, F(A)/F C ^3 (Fig. [9]), one would predict the level of 
polarization to drop to ~0.25 P c at the line center. The relatively 
high signal-to-noise ratio achieved at the Ha line peak (~700 
in the central, 3 pixel bin) allowed us to place an upper limit on 
the depolarization, which must be smaller than 0.2% during both 
epochs. This definitely excludes any depolarization at the level 
predicted by the above relation. 

Finally, we searched for a polarization signal associated with 
the pronounced Ca II H&K absorption, after vectorially subtract- 
ing the adjacent continuum polarization. Although we achieved 
a marginal detection during the first epoch (0.6±0.2%), this sig- 
nal is not present during the second and, therefore, we do not 
attach too great significance to this finding. 



5. Discussion 

The data presented in this paper show that the polarization ob- 
served in NGC 300 OT2008-1 is generated by at least two, al- 
most perpendicular, components Pi and Pi. Because of the con- 
stant polarization angle (rather precisely aligned along the local 
spiral arm of NGC 300), and the wavelength dependence (rem- 
iniscent of a Serkowski law), we have interpreted Pi as being 
caused by ISP. The unusual shape of Pi(X), with its large slope 
parameter K, might be indicative of a dust mixture that differs 
in both chemical composition and grain size from that typical of 
the Milky Way. As of today, there is no consistent model that 
explains the empirical Serkowski law (Draine 2003a), and it is 
therefore impossible to derive the dust properties from the ob- 
served wavelength dependence. 

After reasonably constraining the origin of Pi, what remains 
to be clarified is the nature of P\. 



5.1. Electron scattering 

The most natural source of wavelength-independent linear polar- 
ization in a stellar envelope is Thomson scattering by free elec- 
trons (van den Hulst ll957l ). This is, for instance, the source of 
continuum polarization in supernovae (H6flich [1991| ). More gen- 
erally, polarization by electron scattering is supposed to origi- 
nate in stars with ionized and spatially extended atmospheres, 
such as Be and shell stars, Of and Wolf-Rayet stars, and early- 
type supergiants (Brown & McLean ll977l ). A non-null net polar- 
ization is indicative of an asymmetry in the continuum emitting 
region, and can be understood in terms of partial cancellation of 
the intrinsic polarization generated by the Thomson scattering. 
In the case of NGC 300 OT2008-1, and based on the assump- 
tion of a spheroidal photosphere, the measured level of polar- 
ization (Pi >1.2%) would indicate a substantial asphericity, of 
axial ratio <0.8 (Hoflich 1 199 lb . This geometry would produce 
a wavelength-independent polarization, characterized by a con- 
stant polarization angle (perpendicular to the major axis), similar 
to what is observed. 
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Fig. 8. Best fit solution for 08-07-2008 compared to the ob- 
served Stokes parameters (lower panel), position angle (upper 
panel) and polarization level (upper panel). The curves trace the 
wavelength-independent component Pi (dashed), P2 (dotted- 
dashed) and the global solution (solid). 



Fig. 9. Spectropolarimetry of the Ha line on 01-07-2008. The 
data were binned to 3 pixel (9.6A). The flux spectrum in the 
upper panel is unbinned. The thin colored line traces a high- 
resolution spectrum (AA/A -48,000) obtained with VLT-UVES 
on 02-07-2008 (Pastorello et al. in preparation). 



Evidences of electron scattering polarization (at a level of 
about 0.5%) were found in the early phases of Nova Cygni 
19 92 (Bj orkman et al. [19941 and V838 Mon (Wisniewski et 
al. 120031 1. In both cases this was interpreted in terms of asym- 
metry in the continuum-forming region, identified with a flat- 
tened, spheroidal shell ejected during the outburst. Since in the 
optically thin case the polarization level is proportional to the 
electron-scattering optical depth (Brown & McLean [19771 1. the 
observed decrease in the polarization level is produced by the 
expansion of the shell. The systematic, gray-decline that we ob- 
serve between our two epochs (see Sect. 4. 1 1 might be inter- 
preted along the same lines. 

An important problem with the electron scattering scenario 
in NGC 300 OT2008-1 is the absence of line depolarization 
(Sect. 4.2 1. This effect is expected when the emission lines 



arise in an outer ionized region, where Thomson scattering has 
a lower optical depth (Harrington & Collins 1968; Clarke & 
McLean 1 19741 1. It is observed in a variety of objects, includ- 
ing Be stars (McLean & Clarke IT9791 I. novae (Bjorkman et al. 



[19911 Kaw abata et al. 2006), supernovae (SN 1987A, Cropper 
et al. 1988), eruptive events such as V838 Mon (Wisniewski et 
al. 120031) . and the post-red supergiant IRC+10420 (Patel et al. 
2008 ). However, we note that there are exceptions, in which con- 
tinuum polarization is observed but no line effect is detected, 
both in Be stars (McLean & Clarke 1979; Quirrenbach et al. 
rT997T > and in Nova Sagittarii 1999 (Kawabata et al. 120001 . 

The narrow absorption profile detected on top of the emis- 
sion lines (Bond et al. 2009; Berger et al. 120091 For the case of 
Ha see also Fig. [9] here, upper panel) might enhance the polar- 
ization within the overall line profile. A similar phenomenon has 
been proposed to explain the absence of line depolarization in 
the shell star 77 Cen (McLean & Clarke 1979 ). Given the modest 
depth of the narrow absorption, we however propose this is not 
sufficient to explain the absence of depolarization. Insufficient 
resolution can also be the responsible for the lack of detection of 



line effects. However, we tend to exclude this for our data, since 
the lines are partially resolved. 

Another possibility is that the continuum and the emission 
lines in NGC 300 OT2008-1 arise in similar regions, i.e., that the 
lines do not form very far above the pseudo-photosphere, so that 
both the continuum and the lines are subject to roughly the same 
amount of electron scattering. Although this scenario can only 
be verified through NLTE calculations, we note that this is quite 
unlikely. There is indeed strong evidence that in both SN 2008S 
and in NGC 300 OT2008-1 lines and continuum originate in 
very different regions. While the photospheric temperature and 
radius decrease with time, the emission lines evolve very slowly 
(Berger et al. 120091 Botticella et al. 120091 . This and the very 
high Balmer decrement lead to the conclusion that the emission 
lines are generated by the interaction between the ejected shell 
and the pre-existing circumstellar material, while the continuum 
forms in an inner region (Berger et al. 2009). 

Finally, there is a third possibility, i.e., that the polariza- 
tion source is located outside both the photosphere and the line- 
forming region. This would probably exclude Thomson scatter- 
ing by free electrons as the source of the observed polarization 
and requires an alternative explanation (see next section). 

The absence of electron scattering in NGC 300 OT2008- 
1 can be explained in terms of a spherically symmetric pho- 
tosphere. Although NGC 300 OT2008-1 and IRC+1 0420 have 
similar spectroscopic properties (Smith et al. |2008 ), the post- 
red hypergiant has a continuum polarization that is higher than 
2%, which is indicative of a significative asymmetry in the 
continuum-forming region (Patel et al. 2008). Whatever the rea- 
son for the lack of line-depolarization in NGC 300 OT2008-1, 
it definitely indicates a radical difference between the physical 
conditions in the continuum/line-forming regions of the two ob- 
jects. 
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Fig. 10. Expected polarization wavelength dependence of a 
Milky Way carbonaceous-silicate dust mixture in a single- 
scattering approximation for Ry-5.5 (upper panel) and Ry-3.1 
(lower panel) and two different scattering angles. Dust model is 
from Li & Draine dZOOTT l and Draine (2003b). 



Fig. 11. Minimum polarization solution with an exponential po- 
larization Pi for 2008 July 8. The color-coded segment in the 
lower right panel is the exponential component Pi, while the 
dots trace the corresponding component Pi. The dashed line 
traces the best-fit dominant axis. 



5.2. Dust scattering 

To explain the absence of line effects in Nova Sagittarii 1999, 
Kawabata et al. (2000) proposed that continuum polarization 
originates in a circumstellar dust cloud. The observed wave- 
length dependence (a power law with polarization significantly 
increasing towards the blue) was interpreted as being produced 
by small dust grains. In the case of SN 2008S, Botticella et al. 
(2009) found evidence of an extended dust shell, which is un- 
derstood to be responsible for the large mid-infrared excess ob- 
served in the early phases. This was confirmed by radiative trans- 
fer calculations applied to mid-IR spectra, showing that a sub- 
stantial amount of dust must have survived the outburst (Wesson 
et al. 12009 l l. The existence of a conspicuous dust shell around 
NGC 300 OT20 08-1, p re-dating the outburst, was demonstrated 
by Prieto et al. (2009) based on a mid-IR Spitzer spectrum ob- 
tained three months after the discovery. All of these results infer 
that dust scattering is a perfectly plausible source of continuum 
polarization in NGC 300 OT2008-1. 

One important problem with the application of the dust- 
scattering scenario is the pronounced wavelength dependence 
this is expected to produce. This can be easily understood as 
follows. If there is dust surrounding the object, part of the light 
emitted outside the line of sight would be scattered into the ob- 
server's direction and add to the radiation traveling directly to 
her. Since the scattering by dust produces a high degree of lin- 
ear polarization (see for instance White [T9791 . depending on the 
dust geometrical distribution this can turn into a significant non 
null net polarization. If we consider an infinitesimal dust vol- 
ume element placed at some distance from the central source 
and an incoming unpolarized packet of photons at a given wave- 
length, to a first appro ximatiorljthe fraction of photons scattered 
into the line of sight (polarized perpendicularly to the plane of 
scattering) is proportional to C ext (A) u(A) <E>(/1, 6) (White |1979l . 



where C ext (A) is the dust extinction cross-section, u)(A) is the 
dust albedo, and <f>(A, 6) is the scattering phase function. This in 
turn depends on the scattering angle 6 and the degree of forward 
scattering g{A). Since the albedo is roughly constant in the op- 
tical domain and both C ext and O tend to increase significantly 
in the blue, this implies that the polarization is higher at shorter 
wavelengths. Using the carbonaceous-silicate grain model by Li 
& Draine (2001) and Draine (2003b) and adopting the usual 
parametrization by Henyey & Greenstein ( 119411 ) for the scatter- 
ing phase function, one can estimate this effect more quantita- 
tively. The result is shown in Fig.[l0]for two different scattering 
angles (30 and 70 degrees) and two different Milky Way dust 
mixture^](/?y=3.1, 5.5). As anticipated, the polarization signifi- 
cantly increases at shorter wavelengths, especially in the case of 
lower Ry, and can be fairly well reproduced by an exponential 
law P(A) = P exp|jS(^ - Ao)/Ao]. For /t =4000A, £=-0.7 for 
an average scattering angle of 0-50 degrees (corresponding to 
the typical value of the forward scattering degree in the optical 
domain, i.e., g ^0.6). 

With this result in hand, one can reverse the procedure 
and, using the observed data, deduce the implied component 
Pi. This is shown in Fig. 11 where we have vectorially sub- 
tracted from the data the exponential component P\ (Po=3.7%, 
/?=— 0.7 and 6*1 = 151 degrees). As in the wavelength-independent 
case, we have considered the minimum polarization solution (see 
Sect. |4. \\. As illustrated by Fig. [TT| the resulting Pi is still 

=—0.96): the position an- 



). As illustrated by Fig. 11 
well aligned with a dominant axis (r„ 
gle is 64.7+0.2 degrees and the RMS deviation is 0.31%, i.e., 
slightly higher than for the wavelength-independent solution (cf. 
Table [2]). This is fully consistent with an alignment along the lo- 
cal spiral arm of NGC 300, which is indicative of the interstellar 
nature of this component. 



5 If the dust optical depth is high, then multiple scattering becomes 
important and acts as a depolarizer. See for instance Patat (2005 ). 



6 Wesson et al. {2009 1 have shown that the dust around SN 2008S 
most likely consists of carbon grains. The usage of a carbonaceous- 
silicate mixture here is just illustrative. 
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Fig. 12. Dominant (upper panel) and orthogonal (lower panel) 
polarization on 2008 July 7 after subtracting the exponential 
component Pi. For presentation, the data were binned to 16 pixel 
(57A). The solid curve on the upper panel is a Serkowski law 
with P max =2.6%, /l ma *=4880A and #=5.6. 
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Fig. 13. Best-fit solution for 2008 July 8 compared to the ob- 
served Stokes parameters (lower panel), position angle (upper 
panel), and polarization level (upper panel). The curves indicate 
the exponential component Pi (dashed), Pi (dotted-dashed), and 
the global solution (solid). 



The position angle G\ was fixed by constraining the dominant 
axis to pass through the origin of the Q - U plane (Wang et al 
2003). As done in the wavelength-independent case (Sect. 4.1 



we decomposed P2 along the directions parallel and perpendic- 
ular to the dominant axis, as shown in Fig. [12] for the second 
epoch. The RMS deviation of the orthogonal component P is 
0.15%, which is consistent with the typical measurement er- 
rors and indicates that there are no statistically significant wave- 
length dependencies in the residual polarization. As for Pj, this 
shows a smooth behavior, and can be reasonably well fitted by a 
Serkowski law with P max =4.S%, 4„ fl . t =4250A, and #=3.0. This 
value of K is still large compared to what is typical of the Milky 
Way, but it is significantly smaller than that required by the 
wavelength-independent component (Sect. 



4.1 



The global solution is compared to the data in Fig. 13 



for the second epoch. The fit is reasonably good, even though 
2 (2.5) and the RMS deviation (0.3%) are greater than in 



x 



greater 

Of course, 



4.1 



the wavelength-independent solution (Sect, 
there is absolutely no reason why the wavelength depen- 
dency of dust scattering in the circumstellar environment of 
NGC 300 OT2008-1 should be similar to that of the Milky Way. 
However, this illustrates that the wavelength-independent solu- 
tion, although giving the most accurate reproduction of the data 
with the minimum number of free parameters, is not unique. One 
possible argument against a strong wavelength dependence of Pi 
is that this would imply higher polarization values in both com- 
ponents. In the example that we have just discussed, P2 reaches 
a maximum polarization of about 5% and P\ is ~4% at 4000A. 
These values are quite high, especially considering that the solu- 
tions are undefined up to an additive constant and so the derived 
polarization levels are only lower limits. 

A detailed treatment of the radiation transfer, including light- 
travel effects, multiple scattering, dust properties, and geometri- 
cal distribution (see Patat 2005 ) would be required for a realistic 
modeling. However, a lower limit to the amount of dust required 



to account for the observed level of polarization can be obtained 
as follows. We assume that the dust is placed at a distance r from 
the supposedly unpolarized source and is distributed in such a 
way that all scattered radiation becomes polarized to the max- 
imum possible level P max (attained for G-n/2 and P max ~0.5; 
see White [1979T l. Based on this assumption, if F S (A) is the flux 
scattered into the line of sight and Fq(A) is the flux isotropically 
emitted by the source, then the observed degree of polarization 
is given by 



P(A) 



F + F s 



*(M) oj(A)C ext (A) 

4-nr 1 



Mh 
m H ' 



where Mh is the implied hydrogen mass (for a given gas-to-dust 
mass ratio) and ntn is the proton mass. Substituting the values for 
the fly =3.1 dust mixture (Li & Draine 200 fl Draine l2003bl for 
/l=5500A (w=0.67, g=0.54, C CA ,/i/=5xlO- 22 cm 2 ) and assum- 
ing maximum polarization {G-n/2), one arrives at the expression: 



Mh 



> 3.8 x 10 



" 7 (— ] 
\AUl 



P(5500). 



For instance, if the scattering material is placed at r=1000 AU, 
one needs Mh ~8x10~ 3 M q to reach a 2% polarization. For 
a gas-to-dust mass ratio of 200 (typical of evolved stars with 
carbon-rich dust; see Matsuura et al. [2009 ), this implies a lower 
limit to the dust mass of 4xl0~ 5 M Q . Clearly, if the dust is placed 
at 100 AUp] this lower limit would decrease by a factor of 100. 
As we said, this is only the mass of dust distributed in an asym- 
metric way, because polarization alone cannot constrain the total 
amount of dust. For instance, a spherically symmetric dust shell 
of 10~ 2 M could be present, but this would not produce any 



7 Prieto et al. (20091 set the outer boundary of dust evaporation for 
NGC 3 00 OT2008-1 to be 100 AU, while for SN 2008S Botticella et al. 
(2009) obtained a significantly larger value (~2000 AU). 
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measurable polarization. Prieto et al. (2009) came to the con- 
clusion that about 10~ 4 M of circumstellar dust are required to 
explain the mid-IR excess in NGC 300 O T2008-1. Similar re- 
sults were found by W esson et al. J20091 for SN 2008S, while 
Botticella et al. d20091 > derived a dust mass of ~10~ 3 M Q . Our 
data suggest that a significant fraction of this material (the exact 
amount depends on the dust distance among other parameters) 
must be distributed in an asymmetric configuration. Whether this 
is a torus, a disk, a bipolar outflow (see Bond et al. 2009; Prieto 
et al. 120091) or a clumpy shell superimposed on an otherwise 
spherically symmetric distribution, we are unable to decide on 
the basis of the available data. Spectropolarimetry obtained dur- 
ing earlier epochs might help to differentiate between the various 
possibilities. 

The fiducial date of the outburst is assumed to be April 17, 
2008 (Berger et al. 2009), i.e., 75 days before our first epoch. 
This implies that, at the time of our observations, the expanding 
radiation front is illuminating dust placed at ~10 4 AU perpen- 
dicular to the line of sight. The light curve peak has a width of 
about 100 days (Bond et al. 2009 ), so that the polarized photons 
that we have detected were scattered in a region that extends be- 
tween the dust evaporation boundary and about 10 4 AU. Given 
that at this large distance the incoming flux is a factor of 100 
lower than at 1000 AU, the bulk of the polarization probably 
comes from material located closer to the source, between the 
dust evaporation radius and a few thousand AU. In this scenario, 
the observed decrease in polarization can be understood as be- 
ing caused by the expansion of the radiation flash, coupled to 
the confinement of dust close to the object and possible inhomo- 
geneities in its distribution. 

If the dust is distributed across a sufficiently large distance 
from the outbursting object perpendicular to the line of sight, so 
that the scattering material produces a long integration time, then 
a case B light echo (Patat et al. 2006b ) might become observable 
as the transient fades out. Its evolution with time could help us to 
better characterize the dust distribution. Unfortunately, the lack 
of a scattered light echo detection in the HST images obtained 
more than 4 months after the outburst (Bond et al. 2009) does 
not allow us to place stringent limits on the circumstellar dust 
distance. The HST-ACS spatial resolution (50mas) corresponds 
to -10 s AU at the distanc e of NGC 300 (1.88 Mpc; Gieren et 
al. 2005 ; Rizzi et al. 2006 ). Future HST observations might clar- 
ify whether there is circumstellar dust at distances greater than 
~10 5 AU from the outburst site. Incidentally, the absence of a 
resolved light echo indicates that the transient must be placed at 
a significant distance behind the interstellar dust associated with 
the local spiral arm of NGC 300. 



6. Conclusions 

We have presented VLT-spectropolarimetry of 
NGC 300 OT2008-1 obtained 48 and 55 days after its discovery. 
The continuum polarization is dominated by a component very 
well aligned along the local spiral arm of NGC 300, which 
we have interpreted as originating in interstellar dust. The 
wavelength dependency of the ISP differs substantially from 
what is typical of the Milky Way. Once the ISP is subtracted, 
a substantial signal, reaching about 1.2%, is detected. This 
polarization component is aligned along a completely different 
position angle. 

The simplest explanation of this continuum polarization is 
electron scattering coupled with a significantly asymmetric pho- 
tosphere. Nevertheless, the absence of any line depolarization 
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(very common in sources displaying similar levels of polariza- 
tion) disfavors this scenario and calls for an alternative origin. 

Although we have investigated different possible reasons for 
the lack of depolarization, the most viable physical mechanism is 
scattering by a circumstellar, asymmetric dust cloud. If this is the 
correct picture (also implying that the emitting photosphere does 
not deviate significantly from spherical symmetry), this means 
that a significant amount of dust (>10~ 5 M Q ) must be present 
during the epochs covered by our observations, either because it 
has re-condensed or because it has survived the radiation flash 
produced by the outburst. 

Besides adding to the growing evidence of substantial 
amounts of circumstellar dust in a completely independent way, 
our observations provide direct proof that dust is not only 
present, but that a significant fraction must have an asymmet- 
ric geometry. This implies that one or more asymmetric outflow 
episodes from the progenitor star must have occurred during its 
past history. 

Appendix A: Effects of second order contamination 
on spectropolarimetry 

The purpose of this Appendix is to evaluate the effect of the 
second order on spectropolarimetry when no order-sorting fil- 
ter (hereafter OSF) is used, and to propose a simple technique 
for its correction. 

A. 1. Second-order contamination 

In general, diffraction-based dispersive elements (such as grat- 
ings and grisms) produce replicas of the first-order (FO) spec- 
trum with dispersions that increase proportionally to the or- 
der number. While higher orders can usually be neglected, the 
second-order (SO) spectrum appears as an additive component 
superimposed on the red part of the first order (FO) spectrum, 
with a dispersion that is ~2 times that of the first order. In prac- 
tice, the flux detected at a given wavelength A is the sum of the 
clean spectrum at A and a fraction of the spectrum detected at 
a bluer wavelength, which we indicate as A'. If s(A) is the spec- 
trum one would record on the detector in the absence of SO con- 
tamination, the spectrum which is actually observed s'(A) can be 
modeled as 

s'(A) = s(A) + r 2/1 (A)s(A') (A.l) 

where r2/i(A) is a function that we wish to derive. The first step in 
this procedure is to derive the relationship between A and A'. For 
this purpose, we used archival FORS1 wavelength-calibration 
exposures taken with and without the OSF GG435 (O'Brien 
2008 ). After performing the wavelength calibration and identi- 
fying 6 isolated blue lines and their red replicas, we derived the 
relation between their first order wavelength and their apparent 
wavelength in the second order. The data points are very well 
fitted by a linear relation 

A-a + bA', (A.2) 

where a=-539.5±0.4A and Z?=2.081±0.001. With this result at 
hand, the ratio r 2 /i can be estimated in the following way. We in- 
dicate by S (A) the intrinsic spectrum of an astronomical source, 
by T(A) the transmission function (including atmosphere, instru- 
ment optics and detector efficiency), by F{A) the OSF transmis- 
sion, and by G\(A) and G2(A) the grism FO and SO order effi- 
ciencies,respectively. If s f(A) and s„f(A) are the observed spec- 
tra obtained with and without OSF, respectively, and based on 
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Fig.A.l. Ratio between observed SO and FO as a function of 
wavelength derived for the 300V grism mounted in FORS 1 . The 
upper scale indicates the corresponding FO wavelength A' . The 
solid blue curve is a best fit with a fifth order polynomial. 



the assumption that the order-sorting filter completely eliminates 
the SO (see below), they can be written as 

s f (A) = S(A)T(A)F(A)G l (A) 

s nf {A) = S(A) T(A) G,{A) + S(A') T(A') G 2 {A'). 



Since rj/i, as defined in Eq. A.l is the ratio of the pure SO 
flux at A to the corresponding FO flux at A', this can be written 

as 



s nf (A) - s f (A)/F(A) 

S n f(A') 



(A3) 



The advantages of this formulation are that it is based on mea- 
surable or known quantities and is independent of the detector 
sensitivity function. It can indeed be easily shown that r2/i(A) is 
the ratio of Go to G\ . 

In principle, any bright astronomical target with a contin- 
uum spectrum observed with and without the OSF can be used 
to derive r2/\(A). Nevertheless, because of the lower sensitiv- 
ity of detectors in the blue and since the effect that we wish to 
measure is produced by blue photons, the most suitable choice 
is a hot star. For this reason, and because it is included in the 
FORS1 calibration plan (O'Brien [20081 , we selected the DOp 
spectrophotometric standard starFeige 110 (a^ooo = 23:19:58.4, 
5 2 ooo=-05:09:55.8), which has a ver y blue (B - V=-0. 30) and 
featureless spectrum (Hamuy et al. [19921 [19941 . The FORS1 
archival data that we selected were obtained on 2008 October 
12 with the grism 300V, with and without the OSF GG435, 
which has a cutoff wavelength Af ~4300A. As for the filter 
transmission curve, we used the tabulated data provided by the 
Observatorjj^] One needs to mention that the GG435 OSF does 
not prevent a small portion of the SO being produced in the red 



edge of the spectrum. Its blue cut-off corresponds to a SO wave- 
length of ~8500A, so that redward of this value a contribution 
from photons bluer than 4350A is still present. This makes the 
derivation of r2/\(A) uncertain between 8500A and the red edge 
of the spectral range covered by the detector. 

After wavelength calibration, we extracted the two spectra 
and applied Eq. A. 3 to derive r2/\(A). The result is presented in 
Fig. - 



A. 1 1 This figure shows that the second-order contamination 
can be significant for very blue targets: for instance, the SO con- 
tamination at 6500A is about 35% of the flux detected at 3380A. 
Clearly, the exact amount of contamination depends on the in- 
put spectrum, which is in principle unknown. Nevertheless, one 
must notice that for the SO contamination only the knowledge 
of the blue part of the spectrum is needed, where the observed 
and the uncontaminated spectrum are identical (s'(A) = s(A)). 
The SO begins to contaminate the spectrum at a wavelength A c 
which corresponds to the detector blue cut-off wavel ength A/,. 



A.2 



For FORS1, it is At, ~3200A, which, according to Eq. 
responds to /1 C =6120A. Below this wavelength, the SO can be 
considered negligible and so one can set r2/\(A)-Q for A < A c . 

Therefore, knowing r2/i(A) one can recover the uncontami- 
nated spectrum by simply inverting Eq. A. 1 and using Eq. A.2 



e.g.. 



, I A - a\ 
s{X) = J{X)-r m {A)s'\—^y 



(A.4) 



A way of expressing the severity of the SO contamination is to 
normalize its absolute amount to the flux at any given wave- 
length, by introducing the quantity A(A) 



A(A) 



s'(A) - s(A) 
s'(A) 



»*2/lW) 



s'(A') 
s'(A) ' 



(A.5) 



which is the fraction of the observed flux produced by the SO. 

The application of this method to the data discussed in this 
paper is illustrated in Fig. A.2 To avoid introducing additional 



noise, we replaced the computed r2/i(A) with a fifth order poly- 
nomial fit to the data (solid curve in Fig. \A.l\ . For our object, 
the contamination reaches a maximum relative value of about 
A=5% at about 8000A which, for most purposes is negligible, 
being smaller than the typical accuracy one can achieve in the 
flux calibration. 

Nevertheless, since polarimetry is in most cases dealing with 
flux differences of a few percent, the impact of SO contamination 
needs to be evaluated. 

A.2. Effects on spectropolarimetry 

In dual -beam instruments such as FORS1, spectropolarimetry 
is performed by evaluating the so-called normalized flux ratios 
Fj(A), derived from the ordinary and extraordinary beams pro- 
duced in the same frame by the Wollaston prism (see Patat & 
Romaniello 2006a for a general introduction). For the canonical 
HWP position angles 6>, = |z (where z'=0, 1, 2, 3) and for an 
ideal polarimeter, these ratios are directly related to the Stokes 
parameters: Fq — Q,F\ — U, F2 = —Q, and Ft, = —U. Since the 
observed fluxes in the ordinary {/!' ) and extraordinary (/)'' £ ) 
beams can be modeled with Eq. 
ratios 



A.l the observed normalized 



http://www.eso.org/sci/facilities/paranal/instruments/fors/inst/Filters/ 



f' i '°(A)-f i ' E (A) 
f!'°(A) + f;- E (A) 
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Fig. A.2. Lower panel: comparison between the observed spec- 
trum of the OT in NGC 300 (dotted line) and the corrected spec- 
trum (solid line). The dashed line traces the r2/\{A) function. 
Upper panel: the relative SO contamination A(A). The results 
above 8500A are uncertain (see text). 



after some manipulation can be expressed as 



F'M) 



Fi(X)+gi(A.) Fj(A') 



(A.6) 



where we have defined 



f°(A) + f, E (A) 



ryiiA) = 



Si(A') 
Si(A) 



rw(A) 



Equation A.6 shows that, in general, the effect of SO con- 



tamination on the final polarization is not obvious, since it de- 
pends on the exact spectral energy distribution and polarization 
of the source. However, when the SO contamination is moder- 
ate (A(A) <0.1), the effect is more readily understood. In those 
circumstances, it is given by s(A')/s(A) * s'(A')/s'(A), so that 
g(A) » A(A) (see Eq. A. 5 1 and one can write 



F'M) 



F i (A) + A(A)F i (A') 
1 + A(A) 



Fi(A) + A(A) Fi(A'). 



In our case, where A reaches a maximum value of ~0.05 and the 
polarization in the blue is of the order of 2%, the absolute effect 
on the derived Stokes parameters is of the order of 0.1%, which 
is comparable to the maximum accuracy one can reach with in- 
struments such as FORS1 (Patat & Romaniello 2006a). For this 
reason, the SO effect on the data presented in this paper can be 
neglected. A more rigorous calculation performed by correcting 
each of the ordinary and extraordinary beams following the pro- 
cedure outlined in Sect. IA.1 Heads to the same conclusions. 



Appendix B: The effect of multiple weak polarizers 

In general, if So(/o, Qo, Uo, Vo) is the Stokes vector that defines 
the polarization status of an incoming beam, its passage through 



an optical system can be described by Si = M ■ So, where Si 
is the emerging Stokes vector and M is the Mueller matrix that 
characterizes the system (Chandrasekhar 1950). To a first ap- 
proximation, a dust cloud can be treated as a partial polarizer. 
The Mueller matrix for a non-birefringent, non-scattering, par- 
tial polarizer with horizontal transmission axis (6=0) is given by 
the expression 



M(0) 




(B.l) 



where k\ and ki (k\ > ki\ 0< k\,ki <1) are the principal trans- 
mittances (see for instance Keller 2000, Eq. 3.24). For &i = l and 
&2=0, one obtains a perfect linear, horizontal polarizer, while 
for k\ = &2=li Eq. B.l implies an isotropic, non-absorbing 



medium. The general expression for a generic position angle 8 
can be obtained by computing the product R(-0) ■ M(0) ■ R(0) 
(Chandrasekhar [19501 ), where R(6) is the usual rotation matrix 



R(0) 



io o o^ 

cos 26 sin 26* 
-sin26» cos26» 
1 



After defining s = sin 26, c = cos 26, k = (k\ +ki)l2, and p=(k\ ■ 
ki)l{k\ + £2), one can write the matrix product as 



M(0) = k 



pc 



( 1 

PC 

ps (1 - VW^s 




c 2 + s 2 



ps 

(1 - yi - P 2 )sc 

s 1 + c 2 ^p 2 










where < k,p < 1. For p=0, one clearly obtains an isotropic 
medium with transmittance k, while for p=l the previous equa- 
tion provides the general expression of an ideal linear polarizer 
with position angle 6 (see for instance Keller 2000; Eq. 3.26). 

Using the above matrix, one can compute the resulting 
Stokes vector to be 

IJk = 7 + Q Q pc + U ps, 

QJk = I oP c + Q (c 2 + s 2 Vl - p 2 ) + U (l - Vi - p 2 ) sc , 
Ui Ik = hps + Qo[ \ - yi - p 2 ) sc + U (s 2 + c 2 yjl-p 2 ) , 

Vi/k = Vo yi - p 2 . 

If the incoming beam is unpolarized, the resulting Stokes vec- 
tor is Si(kIo, klop cos 26, klop sin 20, 0), so that the polarization 
is p and its position angle is 6. If the beam is polarized (either 
because the source is intrinsically polarized or because the beam 
passed already through a polarizing system such as a dust cloud), 
then the situation becomes more complicated. However, when 
the incoming polarization and p are far smaller than 1, the pre- 
vious expressions can be approximated to be 

h ~ kl 

Qi ~ kQ + kl p cos 20 
Ui x kU l) + kI psm26 
Vi * kV Q . 

Therefore, in the case of weak linear polarization, the output nor- 
malized Stokes parameters are simply given by 

< Qi > ~ < Qo > +pcos26, 
<U\> ~ < Uq > +psin26. 
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Fig. B.l. Example composition of two Serkowski laws: 
P max =\.Wc, ,W=3500A, 0=70° (K=0.59) and P max =1.5%, 
/U*=9500A, 61=150° (K=1.59). A Gaussian noise (cr=0.1%) 
was added to the resulting points. The wavelength step is 100A. 

This implies that the effect of multiple weak polarizers (such 
as several dust clouds) can be described as a sum of terms, each 
characterized by a polarization pi and a position angle 6*, . In other 
words, when the polarization is weak, the components can be 
considered to be vectorially additive]^] 

An interesting consequence follows from this conclusion. If 
we express the polarization wavelength dependency as p(A) = 
Pmaxf(X) and we imagine a beam of unpolarized light passing 
through a system of interstellar clouds of various p max and 9 
(with d9/dA-0), then the resulting normalized Stokes parame- 
ters are 

< Q > = 2 Pmaxj fi(X) COS 26/, 

< U > = £ Pmaxj fi(X) Sin 20;. 

Therefore, if the dust in the various clouds produces linear po- 
larization with the same wavelength dependency (i.e., dust prop- 
erties are common to all clouds), then the resulting polarization 
position angle 9 = l/2arctan({//(2) is wavelength-independent 
and the entire system is indistinguishable from a single cloud. 

On the other hand, if the different clouds produce polariza- 
tions that do not obey the same law f{A), then the resulting po- 
larization position angle 9 will be wavelength-dependent. Once 
plotted on the Q, U plane, the polarization data will not lie on a 
straight line passing through the origin and will possibly extend 
over different quadrants. 

An example of a system consisting of two clouds both obey- 
ing the empirical Serkowski law (Serkowski et al. 119751 1. with 
two different values of the polarization peak wavelength A max is 
shown in Fig. B.l 
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